Competition between Kondo and RKKY exchange couplings in Pui ^Am^ alloys 
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To clarify the role of the Kondo effect in screening local magnetic moments of Plutonium 5f- 
electrons as well as its competition to the RKKY interactions we use a combination of density 
functional theory with static Hartree Fock and dynamic Hubbard 1 approximations to calculate 
the strength of both the Kondo exchange, Jk, and of the RKKY exchange, Jrkky, couplings 
for Pui-xAniz system as a function of x. We find that Jk increases despite the atomic volume 
gets larger with the Am doping due to unexpected enhancement of hybridization between / and 
conduction electrons in the vicinity of the Fermi level. At the same time, the RKKY exchange is 
shown to reduce smoothly with increasing x. Our results imply that the Kondo effect should be 
robust against the increase in interatomic spacing of this alloy. 

PACS numbers: 71.27.+a, 71.20.Gj 



Magnetic behavior of metallic Plutonium is a topic of 
hot debates in the current literature 0. Naively one 
expects that Pu f-shcll filled with 5 electrons carries a 
total (spin+orbital) momentum which depending on var- 
ious spin-orbit coupling schemes and crystal field effects 
should result in some non-zero values of magnetic mo- 
ment. This has been confirmed by marry state-of-the- 
art electronic structure calculations @, H, 0, [E] based on 
density functional theory (DFT) in its local density and 
generalized gradient approximations jLDA and GGA) [fj 
and by methods such as LDA+U 0, [1] allowing to incor- 
porate effects of on-site correlation energy U using static 
Hartree-Fock type approximations for the f-electron self- 
energies Q . Experimentally, however, none of the six Pu 
crystallographic allotropes shows local moment forma- 
tion: their spin susceptibilities are small, temperature 
independent and display Pauli-like behavior 10]; spe- 
cific heat measurements indicate absence of magnetic en- 
tropy [fH ; muon experiments did not detect the moment 
with the accuracy of 10 -3 hb 13; combinations of neu- 
tron elastic [l3[ and inelastic |l4j scattering data show no 
convincing evidence in existence of magnetic moments ei- 
ther. 

The apparent discrepancy between theory and exper- 
iment can be attributed to miscounting the number of 
f-electrons and can assume the true configuration to be 
f 6 as it was pointed using a variant of the LDA+U calcu- 
lation [l5| with a different prescription to determine the 
position of the f-band (so called choice of double count- 
ing potential). However, this would produce a completely 
inert f-shell and no large effective masses for the Fermi 
surface electrons which contradicts sharply with specific 
heat data exhibiting an enhanced Sommerfeld coefficient 
[l|. Many-body based calculations [l6| based on a com- 
bination of density functional and dynamical mean field 
[H theories (so called LDA+DMFT [3) have recently 
highlighted another effect: the Kondo coupling of Pu f 
electrons which fluctuate between f 5 or f 6 atomic like 
states with the conduction band made of sd electrons 
producing the Kondo singlet. 



To understand this puzzle, one can try to increase Pu 
atomic volume in order to reduce the effect of hybridiza- 
tion and thus to decrease the value of the Kondo cou- 
pling Jk- As a result, if the Kondo screening mecha- 
nism is in play one expects that at some critical inter- 
atomic distance the local moment would eventually show 
up. Indeed most recent studies[li| of stretched Pluto- 
nium lattice based on the LDA+DMFT calculations with 
Continuous-Time Quantum Monte Carlo (CT-QMC) 
method for the solution of impurity problem [20;] have 
detected that at volumes of the order of 30% larger than 
the volume of the fee 8— phase the temperature behavior 
of spin susceptibilities turns from Pauli-like to Curie-like. 
To simulate similar stretch experimentally, Pui_ I Am I 
alloys have been made, which, however, on the basis of 
magnetic susceptibility, electrical resistivity, and photo- 
electron spectroscopy studies [2ll . I22I [23[ , indicated that 
the character of the 5f states does not vary with the Am 
doping. This came into a contradiction with the DFT 
based study of this system [2~4j. 

We thus see that deducing the behavior of the Kondo 
exchange energy as a function of x is an interesting the- 
oretical problem which may shed a light on the descrip- 
tion of the recent experimental results in the Pu-Am al- 
loy. As a minimal model applicable for this description 
is the model of the Kondo lattice, it is also interesting 
to understand an approximate location of this system in 
the Doniach phase diagram [25| where the competition 
between Kondo and RKKY interactions may lead to ex- 
citing phenomena of quantum criticality and exotic su- 
perconductivity. This, for example, is seen in a different 
class of Pu-based 115 materials such as PuCoGas and 
PuRhGa 5 

In the present work, both the Kondo coupling strength, 
Jk, and the hypothetical magnetic interaction strength, 
Jrkky, are calculated and compared against each other 
as a function of x using supercell based electronic struc- 
ture framework where the f-electrons are treated as lo- 
calized. To deduce Jk we measure the hybridization be- 
tween the / level and the Fermi surface sd electrons. To 
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deduce J rkk y we use a newly developed linear response 
approach [23] based on magnetic force theorem[28] . The 
f-electron self-energies are approximated by their atomic 
Hartree-Fock values using the LDA+U method but the 
results are checked against the LDA+DMFT calculations 
with the self-energies extracted by exact diagonalizing 
many-body atomic Hamiltonians [291 ] . The main finding 
of our work is an unexpected increase of Jk with Am 
doping due to a particular behavior of the hybridization 
function in the vicinity of the Fermi level and simulta- 
neous decrease of Jrkky- The latter is expected as we 
scale up the interatomic distances of the lattice. Our 
calculation shows that Jk always remains larger than 
Jrkky assuming that the Kondo effect is robust against 
the increase in atomic volume of Pui-^Am^ which upon 
doping transforms from the Kondo lattice to the diluted 
impurity limit . It also places this system into the heavy 
fcrmion region of the Doniach phase diagram away from 
the quantum critical behavior. 

Our calculations are performed using the full potential 
linearized-muffin-tin-orbital (LMTO) method including 
relativistic effect of spin-orbit coupling [3(| ■ An effective 
U = 4.5 eV describing the on-site Coulomb repulsion 
among the 5/ electrons is used while the other Slater in- 
tegrals (F( 2 \ and F^) are computed from atomic 
physics, and are subsequently rescaled to 80% of their 
values to account for the effect of screening |3l|. Those 
numbers are known to give a reasonable description of 
the electronic structure for both Pu and Pu-Am sys- 
tem [32|, |33j, [34[. To simulate the effect of alloying in 
our study we used supercells with 4 atoms correspond- 
ing to x = 0, \,\, |, 1 in Pui-^Anij; phase diagram, but 
the results of x > § cannot be taken into account for 
Jrkky because it corresponds to a ferromagnetic (FM) 
phase (x = §) and dilute impurity limit (a; — » 1). We 
have also utilized experimental lattice parameters for var- 
ious x 35]. Our calculated electronic structures are found 
to be consistent with the previous studies [H, HH : the 
5f states of Am are well localized and reside in their 5/ 6 
configuration, corresponding to the filled 5/2 shell; Pu 
5/ states are not fully occupied and found in configu- 
ration with 5.4 electrons; Around the Fermi energy the 
sd conduction electrons prevail through which Pu local 
moments can interact with each other via the RKKY 
mechanism and get screened via the Kondo effect. 

In order to estimate the strength of the Kondo coupling 
we calculate the hybridization function between the 5/ 
and the conduction states, A a p(ijj), which is generally 
expressed via the local Green function for the f-electrons 
as follows [3 

A a p(w) = uo a p - e a/3 - G~p(u>) + S Q/3 (w) (1) 



where 



G Q/3 H = ^( w O k - J ff k -AEM]^ (2) 
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FIG. 1: (Color online) Calculated hyubridization function 
Tr{ImA(u})}/Nd'iT around the Fermi level for Pui-^Ania,. 
Lines correspond to the Am doping with x = 0, j, and i. 



LMTO orbital basis while AS(w) is the self-energy cor- 
rection appeared in the f-electron block only with the 
double counting potential subtracted. The matrix of 
impurity levels e Q( a and the average overlap integrals 
o a p can be found straightforwardly from iJ k ,O k [l~8l ]. 
Here we utilize the static Hartree Fock approximation 
for the self-energy by forcing antiferromagnetically or- 
dered state. This is done to keep the same level of ac- 
curacy with our subsequent evaluations of the RKKY 
interactions. We benchmarked these calculations against 
the Hubbard 1 approximation [2!| assuming magnetically 
disordered solutions and find similar results. Both meth- 
ods do not assume DMFT self-consistency with respect 
to the hybridization function and are reduced to the self- 
consistent determination of charge densities similar to the 
Kohn-Sham procedure in DFT. The imaginary part of 
the hybridization function taken at zero frequency deter- 
mines the strength_of the Kondo exchange according to 
a simple estimate (3(| 



Jk = 



Tr{ImA(0)} U 
irN d N(0) e f (e f + uy 



(3) 



Here H^, Ok, are single particle LDA Hamiltonian and 
overlap matrix written in a general non-orthogonal 



where Nd is the corresponding degeneracy of the model, 
N(0) is the density of states at the Fermi level, and 
ef = Tr{i}/Nd- Thus we see that all the parameters in 
this expression can be evaluated in our calculation where 
we find that the average position of the impurity level 
tf ~l eV and the total density of states iV(0) » 1.5 
states /[eV-atom] at the Fermi level are weakly depen- 
dent functions of Am concentration, and the trend in Jk 
is mainly determined by the behavior of ImA (0). 

Figure Q] illustrates our calculated behavior of 
Tr{ImA(ui)}/NdTr for frequencies around the Fermi level 
and for doping levels x < i assuming full degeneracy 
Nd = 14. Three various lines correspond to the dop- 
ings with x = 0,4, and |. The following conclusions can 
be derived. First, all curves look very similar although 
there is an almost rigid shift of the order of 0.4 eV which 
separates the calculated JmA(u) for various x. Second, 
we monitor the overall trend of decreasing the hybridiza- 
tion between the / and conduction electrons with the 
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Am doping as one can trace the maximum of each plot. 
This is easily understood since interatomic distances get 
larger. However, importantly that for all doping ratios, 
JmA(w) has a valley and a hilltop at the right hand side 
of the valley. For pure Pu, the Fermi level is located close 
to the dip and it gradually climbs up to the hilltop as x 
increases. This results in an unexpected increase in hy- 
bridization for u = 0: Tr{ImA(0)}/N d Tr =0.05 eV for 
x = 0; but jumps to 0.10 eV for x — \. (see Table Q] 
for full compilation of the data). Accordingly, Jk also 
increases because iV(0) remains approximately the same 
for all x values. If one sets Nj, = 14 in Eq. then for 
pure Pu, Jk is 340 K but becomes 660 K for x = j and 
further raises to 770K for x — \ . Then Jk decreases but 
is still sufficiently large as x approaches the dilute impu- 
rity limit. One thus concludes that the Kondo screening 
is the robust effect upon the Am doping which would 
prevent the Pu moment to appear at all x . 

As the description in terms of the Kondo lattice Hamil- 
tonian may be relevant for Pui^^Am^ system, its prop- 
erties should be controlled by the competitions between 
the Kondo and RKKY exchange interactions which, ac- 
cording to the Doniach phase diagram, depending on the 
precise value of Jk may lead to either weakly coupled 
magnetically ordered local moment state or to the Kondo 
screened heavy fermion state in the strong coupling limit. 
It can even put the system in the vicinity of quantum crit- 
ical point where exotic superconductivity is believed to 
occur. In the approximation when only a single conduc- 
tion band hybridizes with the f-level, Jrkky scales sim- 
ply as J^A^O) as seen by using the second-order pertur- 
bation theory for the Coqblin-Schrieffer Hamiltonian[37| . 
It may therefore be expected at first glance that both 
Jrkky and Jk should behave similarly upon doping. 
However, in realistic situations detailed electronic struc- 
ture of the material matters as various intcrband tran- 
sitions contribute to exchange processes and this simple 
trend may be violated. 

In order to estimate the strength of the magnetic inter- 
action between localized 5/ states appeared while map- 
ping the Pu sublattice onto the Heisenberg (pseudo)spin 
Hamiltonian, H = Y] Jf^S^-S^we utilize the magnetic 
force theorem within a rigid spin perturbation method 
[28| . In this framework, Jrr> is given as a second-order 
derivative of the total energy induced by the rotations of 
magnetic moments at sites R and R which can be found 
by calculating the following spin-susceptibility-type in- 
tegral: 



RR d(j> aR d<l)/3 R ' 

= E E — /k+q/ ^ x B ^ k + u"> 

x (k + q/|[axB J i,] /3 |kjV q - (R - R '\ (4) 

Here f^j, <r, and are the Fermi function, Pauli spin 
matrix, and the effective magnetic field at atom _R, re- 



TABLE I: The calculated Tr{7mA(0)}/147r, J K and Jrkky 
for various dopings of Pui-^Ania; alloy. 

Am ratio Tr{ImA(0)}/UTv (eV) J K (K) Jrkky (K) 
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spectively. The latter is given by the difference in the 
electronic self-energies for spin up and spin down elec- 
trons. Since the LDA+U method is employed to recover 
antiferromagnetically ordered state, those become fre- 
quency independent matrices and the evaluation of inter- 
atomic exchange interactions is straightforward [27|. In 
practical calculations using the supercells the total num- 
ber of nearest neighboring J's are different for different 
doping ratios, and we take average values after calculat- 
ing all possible nearest J's. 

We find that our calculated Jrkky exhibits a trend 
opposite to Jr . It decreases as x increases as it is evident 
from Table |TJ The Jrkky in pure Pu [x = 0) is smaller 
than Jk, 134 K, which is reasonable in the sense that it is 
set by the scale J K N(0) and that from the experimental 
standpoint there is no local moment in Pu due to the 
Kondo screening. However, the behavior is quite different 
from a simple trend that Jrkky ~Jk-^(0) as at a; = |, 
Jrkky =100 K, and it becomes 67 K at x = \. This 
must be due to interband transitions presented in Eq.(U]). 
The positive sign of Jrkky refers to the AFM order by 
the convention in Eq. |4]) and it is consistent with the 
assumed AFM ground state. 

The comparison of these two quantities, Jk and 
J rkky i provides us with a clear picture for the magnetic 
properties of Pui-xAmj. It follows that Jk is always 
larger than Jrkky up to x — | which covers up the 
whole range of the experiments up to now [2l|, [12, HH . 
So, if the Kondo screening works for Pu, it should also 
work for the alloy. Moreover, the trend is quite sugges- 
tive as we approach the dilute limit. While we cannot 
extract the value of Jrkky for x > | (for x = |there is 
only 1 Pu atom left in our supercell producing FM solu- 
tion), the overall trend for Jrkky to decrease is expected 
as the inter-Pu distances increase. It is therefore clear 
that Jrkky would decrease further as x approaches to 
unity. The behavior of Jk for large x is controlled by 
ImA (0) as the Fermi level reaches the vicinity of the 
top point of the hybridization function as seen in Fig. [TJ 
Therefore it is beginning to decrease slightly at values of 
x > 7j. Nevertheless, even in the dilute impurity limit, 
where the conduction bands are essentially made of Am 
sd electrons, our calculated Jk does not drop sharply as 
seen from Table |U All this implies that the Pu-Am sys- 
tem is far from the quantum critical behavior and resides 
in the heavy fermion state. 

There are possible sources of errors in our estimates. 
First, the calculated Jrkky may be overestimated by 
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the static approximation, such as LDA+U. This, in par- 
ticular, was found in the previous studies [13, UK of 
transition-metal oxides where the calculated exchange 
interactions depending on the level of approximation for 
the self-energy can be further reduced by about 10-30 
%. Second, the use of more refined impurity solvers and 
corresponding effects of the DMFT self-consistency will 
change our estimated values of Jk, although this effect 
is not expected to be large due to (i) generally small val- 
ues in the f-electron hybridization function, and (ii), the 
position of the f-level at around 1 eV which is pretty 
far from the Fermi energy preventing the extreme sen- 
sitivity of the Kondo temperature. In fact, most recent 
LDA+DMFT studies of this system based on the CT- 
QMC method have confirmed these conclusions [3^ |. 

In summary, using a combination of density functional 
theory with self-energy corrections for the 5f-electrons 
we performed the estimates of the Kondo and RKKY 



exchange couplings for the whole range of dopings in 
Pux-aAnXj alloy. It was found that Jk and Jrkky ex- 
hibit opposite trends: the Jk increases with x which is 
attributed to the details in the behavior of the hybridiza- 
tion function near the Fermi level while Jrkky is found 
to decrease as interatomic distances get larger with dop- 
ing. Comparing these two values provides a clear picture 
of the robust Kondo effect as the origin of non-magnetic 
behavior reported in recent experiments on this system. 
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